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ABSTRACT      f 


There  were  two  objectives  in  the  experiments  performed  on  the 
barotropic  model.     One  was  to  correct  the  lag  of  the  troughs  and 
ridges  noted  in  the  barotropic  model  by  the  Naval  Weather  Research 
Facility,    and  the  other  was  to  introduce  a  baroclinic  development  term, 

The  approach  to  these  problems  was  one  which  could  be  incorpo- 
rated into  the  Fleet  Numerical  Weather  Facility  model  with  ease. 

Included  herein  are  the  theoretical  development  of  the  numerical 
model,    a  comparison  with  the  FNWF  model  and  charts  illustrative  of 
the  results  with  various  values  of  the  pertinent  parameters  for  24 
and  48  hours. 

It  was  observed  that  the  development  term  causes  selective 
development  of  lows  and  highs;  but  when  overdone,    spurious  short 
waves  are  formed  which,    however,    do  not  basically  change  the  major 
wave  location. 

It  was  also  observed  that  the  function  of  latitude  used  here  to 
correct  the  lag  in  the  horizontal  advection  of  vorticity,    resulted  in 
greatly  improved  movement  of  troughs  and  ridges  from  30  degrees 
to  60  degrees  of  latitude,   but  could  be  further  refined  to  produce  even 
better  results. 
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1.       Introduction 

The  numerical  model  for  the  prognosis  of  the  500-mb  level  used 
by  the  FNWF  is  based  upon  the  barotropic  potential  vorticity 
equation  [3} .     This  model  has  been  carefully  "tuned"  and  engineered 
over  a  period  of  several  years  and  is  found  to  produce  quite  satisfactory 
results.     The  use  of  a  numerical  model  derived  from  a  baroclinic 
vorticity  equation,    however,    should  improve  prognosis  by  permitting 
development  and  decay  of  pressure  systems.     Due  to  time  limitations, 
it  was  decided  to  utilize  a  quasi-baroclinic  model  which  makes  use  of 
of  the  initial  ( 1000-500)-mb  thickness  and  a  simplified  prognosis 
thereof  to  compute  a  Sutcliffe  development  term  in  the  500-mb  prog- 
nostic equation.     No  attempt  is  made  to  verify  the  thickness  forecasts 
nor  the  implied  1000-mb  prognosis. 

The  model  developed  and  tested  here  is  still  essentially  a  one- 
parameter  model  and  constitutes  only  an  attemp  to  provide  for  some 
development  in  the  conventional  barotropic  model,    which  is  the  primary 
shortcoming  of  an  otherwise  very  successful  prediction  scheme  for 
mid- troposphere. 

In  addition  to  the  use  of  a  baroclinic  term,    some  experimentation 
was  performed  on  the  horizontal  vorticity  advection  term  of  the 
equation  by  varying  it  as  a  function  of  latitude.     Various  degrees  of 
smoothing  were  also  applied  to  the  prognosis  and  the  effects  on  the 
appearance  and  verification  are  noted. 
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2.       The  FNWF  model    jj] 

The  barotropic  potential  vorticity  equation  is  taken  as 

*  •»  (2.1) 

where  h  is  the  depth  of  fluid  being  considered.     For  application  to  the 
atmosphere,   this  layer  is  considered  to  be  approximately  the  tropo- 
sphere and  to  be  centered  at  the  500-mb  level.     Expanding  equation  2.  1 
and  multiplying  through  by  h  gives 

The  last  term  is  assumed  to  be  composed  of  three  components; 
namely, 

a)  a  mountain  term, 

b)  a  latitude  term, 

c)  a  temperature  term. 
Equation  2.  2  then  becomes 

(0  (a)  to        l«  10  ,2.3) 

Term  (2)  is  the  Helmholtz  term.     It  is  approximated  by 

h VT  ''n       -^   r^t      n^  (2.4) 

where  Z  is  the  height  of  the  500-mb  surface,    f  is  a  mean  value  for  (f  +J   ), 
h  is  a  mean  value  for  h,    and^uu  is  a  constant  determined  by  experiment- 
ation. 
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Term  (5)  is  a  mountain  term.     The  approximation 

-i(W>V*Vhb^.a4-V.VM^--V-7  *4^  (2.5) 

h 

is  made  where  M  is  the  height  of  the  terrain,   and  "a"  is  an  emperically 
determined  constant.     With  other  considerations  and  simplifications 
this  results  in 

St  "R      '  (2.6) 

where  "b"  is  a  tuning  constant  and  H  is  a  (500- 1000)-mb  thickness. 

The  thickness  is  carried  forward  during  the  forecast  a§  an 
independent  field  with  half  of  the  speed  of  the  500-mb  wind. 

In  addition  to  the  above,    other  engineering  and  climatological 
modifications  are  made  in  the  FNWF  model  but  will  not  be  discussed 
here. 
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3.       The  baroclinic  model 

The  baroclinic  model  used  here  is  developed  from  the  horizontal 
advective  model    \Z,   p.    400]    in  which  temperature  changes  are  assumed 
to  result  from  horizontal  advection  only.     In  this  model,    it  is  also 
assumed  that  the  vertical  variation  of  the  magnitude  of  the  thermal 
wind  is  the  same  everywhere  and  the  wind  field  may  be  expressed 

V  =  WA<P)VT  (31) 

where      y7    is  the  vertically  integrated  mean  wind  with  respect  to 
pressure  and    V      is  the  vector  difference  between  the  wind  at  two 
selected  isobaric  levels.     With  the  geostrophic  assumption,     \      is 
the  thermal  wind. 

By  the  same  assumption 

/s/+-A(p)/r  (3.2) 

Integrating  the  vorticity  equation 

£T+V-\7(^>f2^  (3.3, 

after  substituting  from  (3.1)  and  (3.2)  gives 
The  left  side  leads  to 


(3.4) 


:>/,  « 


g+V^(*f)*fl*H.tyy 


(3.5) 
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By  the  geostrophic  approximation 


V-  "f-k*72 


(3.6A) 


j*1~v*jp  <3-6B> 

T 

VT»i-|tH7«  <3-6c> 

/r47VH  (3.6D) 

substituting  (3.  6)  into  (3.  5) 

leads  to 

^ii^aL^jj-ayg^d^^Mt  ?to )      (3  8) 

using 

3.8  becomes 

^7l||-fJ(H,0+AiJ(HJ|-7,-H)  (3.9) 

The  right  side  of  (3.  4)  leads  to  a  Helmholtz  term  for  which  the  FNWF 
version  is  utilized.     Introducing  tuning  coefficients  K     and  K     yields 
to  the  final  operational  equation 


(v^sr  ■  -KiXfo*)  -K&I(H,f  7k») 


(3.  10) 
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where  Z  is  the  500-mb  height,  H  is  the  (500- 1000)-mb  thickness,  and 
A  is  the  FNWF  Helmholtz  term.  Assuming  the  thickness  changes  are 
purely  advective  gives 

where  K„  is  0.  5. 
4 


16 


4.       The  numerical  form  of  the  baroclinic  equation 

Equations  (3.  10)  and  (3.  11)  are  solved  over  a  field  composed 
of  63x63  grid  points  which  are  a  distance  d  apart.     The  appropriate 
map  factor  m  for  the  FNWF  63x63  grid  is  that  for  a  polar  stereo- 
graphic  projection    £lj    •     In  order  to  utilize  the  FNWF  subroutine 
library  and  reduce  computation  time,    fixed  point  arithmetic  was  used. 
The  constants  and  corresponding  scaling  factors  are: 

d  =  381KM  =  3.  81  x  10?  cm 
m  =  2m 

f  =   1.03125  x  10"4  sec"1 

-2 
g  =  980  cm  sec 

D  =  Dx  217    (D  =  Z-Zstd) 


H  : 

A 

=  H  x 

21? 

P     : 

=*« 

2"9 

*    = 

/T  : 

,2"9 

(f?M. 

<$) 

At 

=  1  h 

our  = 

3600secs 

which  result  in 

(7**)*$*-  k;*  ;lk  ,M*k  t($,9)-  K3x^<  .  W*iYh,)$  )         (4.  i) 

which  is  the  form  of  equation  (3.  10)  used  in  the  FNWF  Helmholtz 
subroutine     [lj    . 

The  thickness  is  carried  forward  during  the  forecast  as  an 
independent  field  with  half  of  the  speed  of  the  500-mb  wind 

17 


When  expressed  in  finite  difference  form,    (4.  Z)  becomes 

AHs-;?*.?^*^**  j(bjti)  (4.3) 

Both  (4.  1)  and  (4.  3)  are  computed  on  an  hourly  basis.     The  first 

time  step  is  produced  by  the  method  of  forward  time  differences 

^m+ 1       ^m       j.     rn  ,  A     A% 

D  =  D      +  LD  (4. 4) 

If  continuedjthis  method  would  lead  to  computational  instability.     For 
this  reason,    all  other  time  steps  are  produced  by  the  method  of 
centered  time  differencing  4 

D  =  D  +  2AD  (4.  5) 

In  the  event  a  smoothing  operation  is  performed  on  a  field,    a  forward 
time  difference  is  used  as  the  subsequent  time  step. 
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5.       Research  procedures 

Comparing  equations  (2.  5)  and  (3.  10)  shows  that  the  FNWF  model 
without  certain  refinements  closely  approximates  the  baroclinic  model 
when  K     is  zero.     For  this  reason,    this  baroclinic  model  without  the 
development  term  can  be  referred  to  as  the  quasi-FNWF  model; 
however,   the  FNWF  model  has  additional  refinements  not  included 
here,    hence  comparisons  of  verification  scores  are  not  strictly  valid. 

The  data  used  for  the  experiments  were  recent  500-mb  and  (500- 
1000)-mb  thickness  fields  obtained  from  FNWF.     The  actual  computing 
was  done  on  the  Postgraduate  School  Computer  Facility  CDC  Model 
1604. 

The  results  of  each  experimental  run  were  tested  by  the  pillow 
and  RMSE    £lj    ,    then  analyzed  and  visually  compared  with  the  analysis 
at  verification  time  as  well  as  with  other  experimental  runs. 

il»»l  III     — 

X 


Pillow  _.N»i 


£    Cfa-O-ptu^Li. 


RMSE 


It  is  apparent  that  the  RMSE,    although  a  rapid  method  of 
quantitative  comparison  of  fields,    will  not  show  a  variation  of  develop- 
ment of  systems.     It  was  therefore  necessary  to  use  the  visual  method, 

The  experimental  procedure  consisted  of  varying  the  values  of  K 

or  K     or  both.     The  values  of  K     varied  from  0  to  2.  0  and  also 
3  3 

included  variations  of  the  tuning  constants  with  respect  to  latitude. 
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In  addition  to  the  above,    the  fields  were  compared  after  different 
degrees  of  smoothing  were  applied. 
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6.       Experimental  results 

Consider  first  the  advection  term  only.     The  speed  of  movement 
of  the  troughs  and  ridges  were  modified  by  the  use  of  K     as  a  function 
of  latitude.     The  need  for  this  modification  was  made  evident  by  a 
study  of  the  FNWF  model  conducted  by  the  Naval  Weather  Research 
Facility,    Norfolk,    Virginia,    which  indicated  consistently  slow  move- 
ment of  troughs  and  ridges  through  certain  latitudes,    a  deficiency  that 
appears  to  characterize  barotropic  models  in  general  and  even  some 
baroclinic  models. 

After  several  trials,    an  emperical  correction  was  made  for  this 
lag  by  the  employment  of  a  function  of  latitude  (SF)  as  the  factor 

.  r  +  $  i  *  <f> 

.5*  +-  SH1^ 

applied  to  the  horizontal  vorticity  advection, 

The  graph  of  this  function,    Fig.    1,    shows  a  maximum  of  1.  37 
at  21.  4  degrees  and  a  value  of  one  at  the  equator  and  pole. 

The  effect  of  the  use  of  this  sine  function  can  be  seen  by  comparing 
Fig.    2,    the  24-hour  prognostic  chart  without  the  sine  function  and 
Fig.    3,    in  which  the  sine  function  is  applied. 

There  is  an  obvious  displacement  of  all  systems  while  the 
contour  shape  remains  essentially  unchanged. 

To  test  the  effectiveness  of  this  function  on  independent  data,    a 

test  was  run  on  seven  consecutive  24-hour  periods  from  0000Z  24  March 

to  0000Z  27  March.     The  displacement  of  troughs  and  ridges  was 

measured  in  degrees  of  latitude  and  an  average  deviation  from  the 
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Fig.   3 
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TABLE  I 

24-Hr  Prog 

AVERAGE  DISPLACEMENT    FROM 
ANALYSIS 
IN  DEGREES  OF  LATITUDE 

LAT  WITH  SINE  FUNCTION  WITHOUT  SINE  FUNCTION 

TROUGHS    RIDGES    COMBINED      TROUGHS    RIDGES    COMBINED 

10  -2.72(11)  -2.0(2)  -2.62(13)  -2.64(11)  -2.0(2)  -2.54(13) 

20  1.43(23)         -.36(14)  .76(37)  .61(23)  -.93(14)  .03(37) 

30  .70(36)  -1.91(23)  -.32(59)  -.97(36)  -2.26(23)  -1.47(59) 

40  .51(35)  -2.40(30)  -.83(65)  -.66(35)  -3.57(3CX)  -2.00(65) 

50  -.35(31)  -2.15(26)  -1.16(57)  -1.52(31)  -3.46(26)  -2.40(57) 

60  -1.66(18)         -.95(20)  -1.30(38)  -2.38(18)  -2.15(20)  -2.26(38) 

(#  of  samples) 
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TABLE  II 
24-Hr  Prog 
TEST  OF  SINE  FACTOR 
PILLOW^      RMSE* 


V.  T. 

NO 

SF 

WITH  SF 

FNWF 

00Z  24  Mar 

-1 

93 

-2 

94 

1.  1 

104 

12Z  24  Mar 

-5 

101 

-5 

102 

6.3 

115.  3 

00Z  25  Mar 

-5 

90 

-5 

91 

6.6 

108.  3 

12Z  25  Mar 

5 

97 

5 

97 

-3.  8 

105.9 

00Z  26  Mar 

0 

113 

0 

113 

-    .  3 

111.  1 

12Z  26  Mar 

-2 

117 

-2 

116 

-1.  2 

115.0 

00Z  27  Mar 

18 

94 

18 

95 

-25.  5 

101.  5 
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corresponding  analysis  was  taken  for  the  entire  period. 

The  results  of  this  test  are  summarized  in  Table  I.     The  first 
item  observed  by  this  summary  was  that  the  movement  of  troughs 
and  ridges  was  truly  in  need  of  an  increase  at  almost  all  latitudes 
tabulated.     It  is  seen  that  the  use  of  the  sine  function  produces  a 
distinct  improvement  from  30  to  60  degrees.     This  is  particularly 
evident  in  the  "combined"  column  which  includes  troughs  and  ridges. 

Table  I  also  shows  a  possibility  of  too  much  correction  at  20 
degrees. 

The  results  at  10  degrees  are  considered  to  be  inconclusive  due 
to  the  small  number  of  samples  and  the  difficulty  of  defining  troughs 
at  this  latitude. 

Table  II  compares  the  pillow  and  RMSE  for  the   verification  dates 
of  the  test  without  the  sine  function,   with  the  sine  function  and  for 
the  FNWF  product.     The  use  of  the  sine  function  produces  little  change 
in  this  statistic. 

Turning  now  to  the  development  of  systems,    the  problem  was  to 
define  a  proper  K     term  or  to  determine  how  much  development  term 
to  use. 

The  experimental  data  was  selected  at  random  from  the  FNWF 
files  of  height  fields  and  corresponding  thickness  fields. 

The  24-hour  prognostic  tests  were  run  from  0000Z  25  Jan  1967, 
Fig.    4,    to  0000Z  26  Jan  1967,    Fig  5.     This  date  was  chosen  at  random 
and,    although  not  the  best  of  all  possible  situations  to  illustrate 
development,    does  show  the  effects  of  the  several  variations. 
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Fig.  4 
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Fig.   5 
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For  purposes  of  comparison,    the  following  factors  are  considered: 

a)  the  pillow(P)  and  RMSE(R)  in  feet 

b)  illustrative  charts 

c)  the  mean  of  the  deviation  of  selected  low  and  high  central 
values  in  meters  from  the  analysis,    together  with  the  standard 
deviation  from  this  mean. 

It  will  be  seen  that  a  low  RMSE  does  not  necessarily  indicate  a 
good  verification  and  therefore  must  be  carefully  used  in  conjunction 
with  other  evidence  before  drawing  conclusions. 

For  information,   the  FNWF  pillow  and  RMSE  for  this  date 
were  3.  4  and  115.  9,    respectively.     The  mean  of  the  deviation  of  lows 
was  7.  8  with  a  range  of  92  to  -46  and  a  standard  deviation  of  47.  1, 
and  the  mean  of  the  deviation  of  the  highs  was   11.7  with  a  range  of 
31  to  -2  and  a  standard  deviation  of  14.0. 

Caution  should  be  used  in  comparing  the  FNWF  barotropic  results 
to  the  barotropic  model  used  here,    where  K     is  0,    as  this  model  only 
approximates  the  FNWF  model.     However,    comparison  of  cases  with 
and  without  the  development  term  should  indicate  what  effect  the 
development  term  would  have  on  the  FNWF  model. 

The  500-mb  analysis  during  the  forecast  period  showed  little 
significant  movement  or  development  with  the  exception  of  the  low 
pressure  system  on  the  northwest  coast  of  the  U.    S.    which  moved 
inland  and  filled,    while  the  low  in  the  Gulf  of  Alaska  deepened  by  90 
meters  and  moved  east  leaving  a  trough  trailing  behind. 
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Pig.  6 
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Case  I 

Fig.    6  is  the  result  of  using  no  development  term,    K     =  0.     The 
general  position  and  intensity  of  pressure  systems  are  good  and 
follow  the  history  intensities.     The  low  on  the  west  coast  of  the  U.    S. 
fails  to  fill  while  the  low  in  the  Gulf  of  Alaska  deepens  by  only  20 
me  ters. 

mean  deviation  and  range  standard  deviation  from  mean 

P       -7         lows     -5.0(51,-38)  29.6 

R     105        highs    19.  2(28,  9)  6.  8 

K     was  next  set  at  1.0.     The  prognosis  became  more  unstable 
with  each  time  step  and  after   13  hours,   the  FNWF  Helmholtz  sub 
routine  £Q  would  not  converge. 

Various  other  values  of  K     were  tried.     It  was  found  that  when  K 
is  0.  5,    the  model  runs  for  18  hours,    K    at  0.  3  for  23  hours  and  K 
equal  to  or  less  than  0.  2,    over  24  hours. 

Since  the  hourly  change  of  the  thickness  field  is  based  solely  on 
horizontal  advection,    it  is  believed  that  this  field  becomes  less 
representative  of  the  actual  field  fairly  soon  and  is  probably  a  major 
source  of  observed  instability.     For  this  reason,    K     was  divided  by 
2  every  six  hours  in  the  following  examples  unless  otherwise 
indicated.     Other  periodic  reductions  of  K     were  tried,    however,   this 
one  gave  the  best  results. 

In  addition,   the  FNWF  subroutine  SAU     jjQ    is  employed  on  the  D 
and  H  fields  every  three  hours  to  suppress  the  generation  of  short 
waves. 
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Fig.  7 
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Fig.  8 


34 


Case  II 

In  Fig.    7,    K     is  0.  5.     Although  the  major  features  remain 
unchanged,    compared  to  Fig.    6,    there  is  some  increase  in  short  wave 
development.     The  mid-latitude  lows  are  from  20  to  40  m  deeper  and 
highs   10  to  30  m  higher. 

mean  deviation  and  range      standard  deviation 
P  -8      lows     -8.3(45,-39)  22.7 

R        115      highs  26.4(50,  8)  14.9 

Case  III 

In  Fig.  8,    K     is  increased  to  1.0.     Further  short-wave  development 
is  noticed  in  addition  to  stronger  development  of  both  highs  and  lows. 

The  development  is  selective,    as  can  be  seen  by  the  low  north  of 
Siberia  remaining  unchanged  while  the  low  in  Idaho  deepens  by  50  m 
and  the  one  in  the  Gulf  of  Alaska  by  80  m. 

mean  deviation  and  range        standard  deviation 
P  -8      lows  -16.  8  (2,  -43)  15.7 

R        128      highs  34.6(62,  8)  22.3 

In  the  above  examples,    it  becomes  evident  that  the  use  of  the 
development  term  produces  an  unacceptable  amount  of  short-wave 
development.     FNWF  has  a  subroutine,    SAG    £1J     ,    which  can  be  used 
to  remove  some  or  all  perturbations.     It  is  of  the  form 

A  =  A  +  k  7  A 

The  problem  remains  as  to  what  size  to  make  k  and  how  many 
times  to  use  the  smoother. 
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Fig.  9 
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Pig.  10 
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The  smoother  was  tried  with  k  =  0.  125  and  used  three  successive 
times  on  the  D  and  H  fields  every  six  hours.     The  results  of  this  heavy 
smoothing  in  addition  to  no  sine  function  and  K     at  0  is  seen  in  Fig.    9. 

Case  IV 

Many  of  the  systems  have  been  smoothed  out  of  existence  and 
others  have  become  unrepresentative  of  the  verifying  analysis.     This 
example  clearly  shows  that  a  low  RMSE  does  not  necessarily  mean  a 
good  prognosis. 

mean  deviation  and  range  standard  deviation 

P  -7      lows  116.0(174,  92)  32.1 

R        105      highs    -29.6(5,-77)  34.2 

Case  V 

Maintaining  the  heavy  smoothing,    and  setting  K     at  1.0,    starts  to 
bring  the  prognosis  pattern,    Fig.    10,    back  to  the  analysis,   but  not 
enough. 

mean  deviation  and  range        standard  deviation 
P  -8  lows   110.  3(142,  84)  22.4 

R        107  highs  -26.  2(3, -60)  32.0 

Case  VI 

In  Fig.    11,    K     is   1.  0  but  is  not  diminished  during  the  24  hour 
period.     As  noted  earlier,    this  would  be  too  unstable  to  be  computed, 
were  it  not  for  the  use  of  a  heavy  smoother. 

Compared  to  Fig.    9,    both  highs  and  lows  show  much  greater 
development,    but  this  is  attained  with  an  excessive  gain  in  short 


waves, 
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Fig.  11 
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Pig.  12 
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Fig.  13 
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mean  deviation  and  range         standard  deviation 
P  -8  lows  70.9(117,  14)  41.2 

R        147  highs  8.4(63, -30)  30.5 

Case  VII 

Continuing  the  search  for  a  compromise  between  development  and 
smoothness,    K     is  again  divided  by  2  every  6  hours  and  the  initial 
value  of  K     is  set  at  2.0,    which  results  in  Fig.    12.     While  the  central 
heights  are  closer  to  the  analysis  than  Fig.    9  or  10,    there  is  still  too 
much  short-wave  development,    although  less  than  Fig.    11. 

mean  deviation  and  range         standard  deviation 
P  -8         lows  85.  5  (117,  55)  23.9 

R        128  highs  6.  2(53, -38)  33.0 

Case  VIII 

In  Fig.    13,    the  smoothing  coefficient  is  reduced  to  1/16  vice  1/8 
while  K     is  1.0.     This  is  seen  to  produce  development  with  acceptable 
smoothness. 

mean  deviation  and  range         standard  deviation 
P  -8         lows  48.  9(69,  28)  14.7 

R       116  highs  4.0  (36, -21)  18.8 

Case  IX 

Continuing  in  the  same  vein,    the  smoothing  coefficient  is  reduced 
to  1/32  while  K     remained  1.0,    Fig.    14.      Further  development  is  noted 
but  so  is  considerable  short-wave  formation.     The  compromise  point 
has  been  passed. 
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Fig.  14 
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Fig.   15 
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Fig.  17 
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mean  deviation  and  range         standard  deviation 
P  -8         lows  .  3(26,  -22)  17.  5 

R        140  highs  27.6(65,  7)  21.6 

The  above  findings  are  now  tested  for  a  48  hour  period  using  a 
verification  analysis  of  1200Z  25  Jan.    1967,    Fig.    15.     The  FNWF 
pillow  and  RMSE  for  this  date  were  -16.0  and  173.  2,    respectively. 
The  mean  deviation,    range  and  standard  deviation  of  the  lows  were 
-81.  0(119,  -230),    119  and  for  the  highs  -24.  3(-10,  -39)  and  11.  8. 

Case  X 

Starting  with  the  heavy  smoother:     k  of  1/8,    3  scans  every  6  hours, 
and  K     at  0,    results  in  Fig.    16.     This  is  an  extremely  smooth  chart 
with  most  baroclinic  features  missing.     A  misleading  low  RMSE  is 
again  encountered. 

mean  deviation  and  range         standard  deviation 
P  8         lows  117.0(165,  81)  102.7 

R        143  highs  -154(1   sample) 

Case  XI 

Continuing  with  the  heavy  smoother  and  changing  K     to  1.0  produces 
some  slight  development,    Fig.    17,    but  still  far  short  of  the  analysis. 

mean  deviation  and  range         standard  deviation 
P  7  lows  101.  3(147,  66)  108.4 

R        142  highs  -154(1   sample) 

Case  XII 

Reducing  the  smoothing  factor  to  1/16  with  K     of  0  improves  the 
chart,    Fig.    18,    over  the  previous  two  results  but  still  not  enough. 
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Fig.  18 
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Fig.  19 
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mean  deviation  and  range         standard  deviations 
P  7  lows  67.0(112,  37)  110.  8 

R        148  highs  -126(1   sample) 

Case  XIII 

With  the  smoothing  coefficient  at  1/16  and  K     at  1.0,    further 
development  is  observed  in  Fig.    19,    but  even  more  would  be  desirable. 

mean  deviation  and  range         standard  deviation 
P  7  lows  35.4(84,  -5)  119 

R        161  highs  -67. 0(-24,  -110)  43 

An  attempt  was  made  at  reducing  k  to  1/32  while  leaving  K     at  1.0. 
This  variation  produced  a  30  hour  prognosis  before  overflow  occurred. 

Case  XIV 

The  final  variation,    Fig.    20,    was  the  most  satisfactory  although 
not  spectacular.     In  this  variation,    k  is   1/16  but  used  twice  instead  of 
three  times  every  six  hours  while  K     is   1.0. 

mean  deviation  and  range         standard  deviation 
P  8         lows  -55.4(-14, -123)  136.5 

R        168  highs  5.0(62, -62)  52.4 

The  results  of  the  above  cases  are  summarized  in  Table  III. 
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7.       Independent  data  results 

The  last  combination  of  factors,    that  is,    K     of  1.0  reduced  by  one 
half  every  six  hours  and  a  smoothing  coefficient  of  1/16  applied  to  the 
height  and  thickness  fields  every  six  hours,    was  applied  to  independent 
data  of  seven  consecutive  24-hour  prognostic  charts. 

These  charts  covered  the  period  0000Z  24  March  to  0000Z  27  March. 
A  comparison  was  made  with  full  development  term  of  1.0,    no  develop- 
ment term  and  the  FNWF  product  for  these  dates. 

Upon  summarizing  the  results  of  the  above  tests,    it  was  seen  that 
a  K     of  1 .  0  was  probably  too  high  and  another  group  of  seven  charts 
was  produced  with  a  K     of  0.  7. 

The  results  of  all  of  the  above  are  listed  in  Table  IV.     The  pillow 
and  RMSE  for  the  variations  on  the  applicable  dates  is  presented  in 
Table  V. 

Table  IV  indicates  a  sharp  reduction  in  the  mean  deviation  of  the 
lows  from  the  analysis,    however,   the  choice  of  a  K     of  0.  7  was  too 
much  of  a  reduction  from  1 .  0  in  the  case  of  the  lows  but  was  quite 
good  for  the  highs.     A  K     of  about  0.  8  would  probably  have  produced 
better  results  for  this  period. 

Compared  to  FNWF,  a  K  of  0.  7  results  in  a  good  reduction  in  the 
range  of  deviations  for  both  highs  and  lows  and  a  marked  improvement 
for  the  standard  deviation  from  the  mean  deviation  for  highs. 
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TABLE  V 

INDEPENDENT  DATA  TEST  OF  DEVELOPMENT  TERM 
PILLOW  AND  RMSE  (ft.  ) 

V.  T.  K3  =  0  K3=1.0  K3=.7  FNWF 

0000Z  24  Mar 

1200Z  24  Mar 

0000Z  25  Mar 

1200Z  25  Mar 

0000Z  26  Mar 

1200Z  26  Mar 

0000Z  27  Mar 


p  -1 

-2 

-2 

1.  1 

R  93 

96 

92 

104 

P  -5 

-5 

-5 

6.3 

R  101 

106 

100 

115.  3 

P  -5 

-5 

-5 

6.6 

R  90 

101 

92 

108.  3 

P  5 

5 

5 

-3.  8 

R  97 

101 

95 

105.9 

P  0 

0 

0 

-.3 

R  113 

113 

109 

111.  1 

P  -2 

-2 

-2 

-1.2 

R  117 

119 

115 

115.0 

P  18 

17 

17 

-25.  5 

R  94 

101 

95 

101.  5 
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8.       Conclusions 

The  modification  of  the  horizontal  advection  term  by  a  function  of 
latitude  as  used  herein,    generally  improved  the  movement  of  troughs 
and  ridges  by  reducing  the  lag  between  prognosis  and  analysis. 

It  appears  that  a  refinement  of  this  function  with  a  displacement 
of  the  maximum  from  about  20  degrees  to  about  45  degrees  would  be 
more  rewarding. 

While  compiling  data  for  this  study,   there  appeared  to  be  a 
relation  between  the  zonal  index  and  the  lag  of  the  troughs  and  ridges. 
It  is  recommended  that  further  study  of  this  phenomena  be  conducted 
to  arrive  at  a  more  exact  relationship  and  how  it  could  be  applied  to 
proper  movement  of  troughs  and  ridges. 

It  is  readily    seen  that,    as  the  name  implies,    the  development  term 
produces  development  in  the  form  of  deeper  lows,    higher  highs  and 
also,    unfortunately,    unwanted  short  waves. 

It  results  in  distinct  and  selective  changes  in  intensities  of  systems 
without  materially  changing  the  location  of  these  systems. 

It  showed  marked  improvement  over  the  FNWF  highs  and  a  lower- 
ing of  the  range  of  errors  in  the  lows. 

The  inclusion  of  this  term  to  some  degree  in  an  operational 
barotropic  model  would  be  beneficial. 

It  is  suggested  that  the  most  desirable  approach  would  be  to 
incorporate  these  empirical  modifications  in  the  FNWF  500-mb 
prediction  equation  while  simultaneously  carrying  out  the  FNWF  1000-mb 
prediction  hour  by  hour.     The  resulting  H  field    could  then  be  incorpor- 
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ated  into  the   500-mb  forecasts  each  hour,    which  should  be  far  more 
reliable  than  that  produced  by  simple  advection  of  thickness. 

A  development  term  so  produced  could  be  made  smaller  initially, 
say  K     =  0.  25,    but  would  not  be  reduced  periodically  since  the  thick- 
ness field  would  presumably  stay  in  phase  better  with  the  500-mb 
height  field,  i 

This  scheme  should  produce  improved  long  range   500-mb  and 
surface  prognostic  charts. 
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